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ABSTRACT 

We present the Snapshot Hubble U-band Cluster Survey (SHUCS), a project aimed at char- 
acterizing the star cluster populations of ten nearby galaxies (d < 23 Mpc, half within « 12 Mpc) 
through new F336W (U band equivalent) imaging from Wide Field Camera 3, and archival BVI- 
equivalent data with the Hubble Space Telescope. Completing the UBVI baseline reduces the 
age-extinction degeneracy of optical colours, thus enabling the measurement of reliable ages and 
masses for the thousands of clusters covered by our survey. The sample consists chiefly of face-on 
spiral galaxies at low inclination, in various degrees of isolation (isolated, in group, merging), and 
includes two AGN hosts. This first paper outlines the survey itself, the observational datasets, 
the analysis methods, and presents a proof-of-concept study of the large-scale properties and star 
cluster population of NGC 4041, a massive SAbc galaxy at a distance of w 23 Mpc, and part 
of a small grouping of six giant members. We resolve two structural components with distinct 
stellar populations, a morphology more akin to merging and interacting systems. We also find 
strong evidence of a truncated, Schechter-type mass function, and a similarly segmented lumi- 
nosity function. These results indicate that binning must erase much of the substructure present 
in the mass and luminosity functions, and might account for the conflicting reports on the in- 
trinsic shape of these functions in the literature. We also note a tidal feature in the outskirts 
of the galaxy in GALEX UV imaging, and follow it up with a comprehensive multi-wavelength 
study of NGC 4041 and its parent group. We deduce a minor merger as a likely cause of its 
segmented structure and the observed pattern of a radially decreasing star formation rate. We 
propose that combining the study of star cluster populations with broad-band metrics is not only 
advantageous, but often easily achievable thorough archival datasets. 

Subject headings: surveys: SHUCS — galaxies: individual (NGC 4041) — galaxies: star clusters: general 
— galaxies: interactions — galaxies: star formation — galaxies: groups: individual: LGG 266 



^Australian Astronomical Observatory, PO Box 915, 
North Ryde NSW 1670, Australia; iraklis@aao.gov.au. 

^Department of Astronomy &; Astrophysics, The Penn- 
sylvania State University, University Park, PA 16802, USA. 

^ Space Telescope Science Institute and European Space 
Agency, 3700 San Martin Drive, Baltimore, MD 21218, 
USA. 

*Max-Planck-Institut for Astronomy, Konigstuhl 17, D- 
69117 Heidelberg, Germany. 

^Departement de Physique, de Genie Physique et 
d'Optique, and Centre de Recherche en Astrophysique du 



Quebec (CRAQ), Universite Laval, Quebec, Canada. 

® Department of Astronomy, University of Wisconsin- 
Madison, 5534 Sterling, 475 North Charter Street, Madison 
WI 53706, USA. 

''Excellence Cluster Universe, Boltzmann-Strasse 2, 
85748 Garching bei Miinchen, Germany. 

^Astrophysics Research Institute, Liverpool John 
Moores University, Egerton Wharf, Birkenhead, CH41 
ILD, UK 

^European Southern Observatory, Karl-Schwarzschild- 
Strasse 2, 85748 Garching bei Miinchen, Germany. 

^"^Department of Astronomy, Stockholm University, Os- 



1 



1. Introduction 

The launch of the Hubble Space Telescope 
{HST) over two decades ago started a revolu- 
tion in the study of extragalactic star clusters. 
The discovery of large numbers of young com- 
pact clusters in star-forming galaxies led to the 
suggestion that they could be the present-day 
analogues of gl o bular clusters (see reviews of 
Whitmor^ 2003 : Larsen 2004a ). The question 
of whether these young clusters can survive for a 
Hubble time has still not been settled but their 
longevity appears to be critically dependent on en- 
vironmental conditions within their host galaxie s 



(|de Griis fc Parmentiej2007tlBastian et al.l201l[) 

The instahation of Wide Field Camera 3 (WFC3) 
on HST has vastly upgraded the imaging capa- 
bilities of the telescope shortward of 4000 A. It is 
now much easier to measure the age and mass dis- 
tributions of large populations of star clusters in 
galaxies, and address fundamental questions such 
as their long term survival chances. To obtain 
ages and extinctions for clusters younger than 
~ 2 Gyr, it is essential t o obtain photometr y 
across the UBVI baseline ( Anders et al ] |2004al) . 
Prior to WFC3, HST U band imaging of suf- 
ficient depth and spatial coverage was feasible 
for only a f ew regio ns of nea rby galaxies (e. g. 
Smith et al.i r2007: Anders et al.. 2004b). and only 



distant systems in their en tirety (e. g. lOstlin et al 
20031: lAdamo et al.|[2010al) 



All local, late-type giant galaxies host popu- 
lations of young and intermediate-age star clus- 
ters, often with masses and densities that ri- 
val globular clusters. It has been proposed that 
the vast majority of stars are formed in clusters 
but that most of these clusters (~ 90%) rapidly 
dissolve (Lada & Lada 2003) due to the expul- 
sion of residual gas from star formation. Clus- 
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ters that survive this period of "infant mortal- 
ity" can disrupt through stellar evolution, two- 
body relaxation and the tidal field of the host 
galax y on Gyr timescales ([Bastian fc Gieles 
l2008h . There has been much debate in the 
literature over whether infant mortality exists, 
whether the later phases are mass-dependent 
or not (e.g. iLamers et al.ll2005t iFall et al.ll2005t 



Whitmore et al.ll2007l : IChandar et al.ll2010l : lBastian et al 



201 II ). a nd even if most star s do indeed form in 
clusters ( Bressert et al. 2010l) . By conducting a 



survey of local star-forming galaxies, it will be 
possible to obtain large samples of clusters cov- 
ering a wide range of masses and ages, and thus 
help answer several open questions. 

In addition, the local environment has recently 
been suggested as a major contributor to clus- 
ter disruption (e. g. lElmegreen fc Hunted l20ld 
Kruijssen et al.l 120111 ). The relationship between 
the formation of clusters and the properties of the 
host galaxy is, in fact, far from clear on both global 
and local scales, and a survey would help to inves- 
tigate these relationships. So far, detailed studies 
of large num bers of clusters in s r nall samples of 



2010: Adamo et al. 2011: Silva- Villa fc Larsen 



galax i es (e. g. iMeurer et al.l 119951: iGoddard et al 



20111 ) or small nu mbers of cluste r s in larg e sam- 
ples of galaxies (LarsenI 2004b; Bastian 2008t 
MuUan et al.l 120111 ) have suggested that environ- 
ments with a higher star-forming density form a 
higher fraction of their stars in clusters. This con- 
nection is also relevant to the link between galaxy 
interactions and increased star formation. Since 
star clusters can trace bursts of star formation, 
they can then provide a viable chronometer for 
the interaction history of a galaxy. 

Naturally, understanding the physics that gov- 
ern these clusters is essential to utilizing them 
as tracers of star formation. After two decades 
of 77S'T-driven research, several cluster parame- 
ters are considered as standard, such as their dis- 
tributions of luminosity and size. Perhaps most 
notably, many studies have investigated the star 
cluster mass distribution, with a common finding 
that it can be represented by a power law of index 
near - 2 Ce. g. IZhang fc Falll Il999l : Ide Griis et aT 



2003: Bik et al 



2003f) . More recent studies have. 



however, suggested that the true form is that of a 
Schechter function, a power-law distribution with 
an exponential truncation at high masses. This re- 
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Fig. 1. — Colour composite imaging of NGC 4041, using all available HST imaging {UBVI and Ha). 
The imprints of various detectors are visible here, arising from uneven spatial coverage across the optical 
baseline. The new F336W imaging envelops the archival WFPC2 pointings, while the ACS image leaves a 
trace of its chip gap. The face of NGC 4041 exhibits a composite structure. Pink Ha bubbles trace the 
distribution of young clusters along spiral arms that wind tighter in the inner regions than the outer galaxy. 
A bright central component seems to define this two-step structure, while the central peak is offset from the 
geometrical center, the center of the outermost isophote, by w 1". 



fleets the mass function of Giant Molecular Clouds 
(e. g. ISolomon et al.lll987| ). out of which star clus- 



ters form. The slope also appears to be a func- 
tion of brightness, that is to say, the brigh ter the 
subsample, the steeper the extracted slope (|Gieles 
20101) . In addition, preliminary indications sug- 
gest that the characteristic 'Schechter mass', M*, 
where th e truncation occurs, depends on envi- 
ronm ent ( Larsen 120091: iGieledboOOt iBastian et al 



20121 ). Again, a large survey of clusters in a variety 
of host galaxies will be able to settle this issue. 



Another parameter of interest is the star cluster 
size distribution and whether this is related to en- 
vironment and/or age. Observations indicate that 
effective radii are constrained to a range of 0.5 - 
10 pc despite the large dynamic range in cluster 
mass (IPort egies Zwart et al. 2010l ). The observed 
radius distribution is well described by a log- 
normal distribu tion with a peak at 3-4 pc for both 
young clusters ( Barmbv et al.|[2"006t Bastian et al 



i2011) and old globular clusters (e. g. iJordan et al 
120051 ) ■ with the exception of the "faint- fuzzy'' 
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clusters discovered by iLarsen fc Brodie (I2OOOI) . 
There are indications that the cluster core ra- 



dius increases with age (e. { 


I. Mackev & Gilmore 


2003aibl: Scheeomaker et al. 


20071 Bastian et al. 


2008;lHurlev & Mackev 2010l). A survey of cluster 



sizes in different environments within galaxies will 
permit studies of the cluster size distribution as a 
function of age and environment. 

From the above, it is clear that an extensive 
survey of a large sample of star clusters in a va- 
riety of environments, aimed at measuring their 
age, mass, and size distributions, will address 
many fundamental questions regarding their prop- 
erties, survival rates, formation histories, and re- 
lated environmental dependencies. In this paper, 
we describe such an endeavor: the Snapshot Hub- 
ble U-band cluster survey (SHUCS) combines new 
WFC3 F336W (U band-equivalent) imaging with 
archival HST _B VT-equivalent imaging, to mea- 
sure the properties of large samples of young clus- 
ters in nearby (mostly c? < 12 Mpc) star-forming 
galaxies. We present a full description of the sur- 
vey in terms of sample definition and data re- 
duction (Section [2]); the detection and photom- 
etry of star clusters, and the derivation of age, 
mass, and extinction (Section [3]). In the second 
part of this paper, we present a proof-of-concept 
study for NGC 4041 a b right (g = 12.1 mag. 



de Vaucouleurs et al.l Il99ll ) face-on SAbc galaxy 
with a double component disk (see Figure [T]). In 
Sections m and [SI we discuss the large scale proper- 
ties and the cluster populations of NGC 4041. We 
summarize this work, as well as our expectations 
for the full survey in Section |6l 

2. The Survey: Target Selection and Data 
Reduction 

To define our galaxy sample for WFC3-F336W 
imaging, we selected galaxies within 25 Mpc that 
have BVI images available in the HST archive. 
We restrict the sample to any galaxy imaged with 
the Advanced Camera for Surveys (ACS) or the 
Wide Field and Planetary Camera 2 (WFPC2), 
but imposed no constraints on image depth. We 
chose this distance limit and instrument set to en- 
sure that individual star clusters are resolved and 
that deep WFC3 F336W imaging can be obtained 
with a 30 min exposure time. The depth is quan- 
tified as reaching mpssew ~ 26 mag at acceptable 



error levels, as will be demonstrated in Section [3?T] 
We excluded dwarf irregular gal axies as they typi - 
cally host very few clusters (e. g. Seth et al. 20041 ). 
but do not impose a strict lower mass limit on the 
dataset. 

The resulting sample of 22 galaxies was drawn 
from the combined target lists of th e 11 Mpc Ha 
Ultra violet Galaxy Survey (llHUGS: lKennicutt et al 



20081): the Loc al Volume Legacy survey (LVL 



Dale et al.ll2009l) : the ACS Nearby Galaxy S urvey 



Treasury (ANGST; iDalcanton et al.l [2009'); the 
HST H a Snap s hot su rvey (HHsnap; JIo 20031); 
and the [LarsenI (|2004bf ) catalog. The last catalog 
on that list contained size measurements, photom- 
etry and information on morphology for clusters 
in galaxies with WFPC2 imaging in various filters 
(typically B, V and I equivalents) and ground- 
bas ed UBV RI imaging. Like the present work, 
the [LarsenI study was also aimed at studying star 
clusters and their immediate environment. 

Our Cycle 18 Snapshot program ''HST U band 
Survey of Star Clusters in Nearby Star-Forming 
Galaxies" (PI Smith, ID 12229) was awarded 22 
snapshot orbits, one for each of 22 proposed tar- 
gets. Eleven galaxies were observed, as listed in 
Table [1] as per the nominal completion rate of 
50% for SNAP programs (one observation failed) . 
About half of observed galaxies are within « 
12 Mpc, and four systems lie at 18 — 23 Mpc. All 
exposures were taken with the F336W filter and 
the duration was 1800 s. A three-point dither line 
pattern was chosen to cover the chip gap and to 
aid in the removal of hot pixels, cosmic rays and 
other artifacts. Either the UVIS-FIX or UVIS2- 
FIX aperture was used depending on the spatial 
extent of the galaxy compared to the 162" x 162" 
UVIS field of view. The precise pointings were 
chosen to give maximum overlap with the archival 
BVI observations. 

Throughout this series we will use the Johnson 
filter notation with the specific HST filter sub- 
scripted, but will not at any point convert between 
the two systems. For example, F336W will be 
denoted as Usse, while F555W and F606W, both 
roughly corresponding to Johnson V band, will be 
written as V555 and Vgoe- The C/336 through /814 
base will be referred to as UBVI. 
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Table 1 
WFC3/UVIS Target List. 



Name 


RA 


DEC 


Morphology 


Distance 


log (A/.) 


Rcf. 


Other HST data 




(h m s) 


(o / //) 




(Mpc) 


(Mo) 






NGC 247 


00 47 10.49 


-20 46 09.00 


SAB(s)d 


3.6 


9.06 


[1] 


FllOW, F160W(N) 


NGC 672 


01 47 54.06 


+27 25 55.80 


SB(s)cd 


8.1 


9.43 


[2] 


F658N(A) 


NGC 891 


02 22 32.90 


+42 20 45.80 


SA(s)b? 


10.2 


10.84 


[1] 


F250W(A), F656N(W2), F160W (N) 


NGC 925 


02 27 05.14 


+33 34 54.50 


SAB(s)d 


9.3 


9.65 


[3] 


F160W(N) 


NGC 1003 


02 39 16.40 


+40 52 20.40 


SA(s)cd 


11.1 


9.44 


[3] 




IC 356 


04 07 46.47 


+69 48 45.20 


SA(s)ab pec 


11.2 


11.32 


[2] 


F658N(A) 


NGC 2146 


06 18 37.71 


+ 78 21 25.30 


SB(s)ab pec 


17.2 


11.04 


[4] 


F658N(A), F160W(N) 


NGC 2997 


09 45 38.70 


-31 11 25.00 


SAB(rs)c 


12.2 


10.61 


[2] 


F220W, F330W(A) 


NGC 3756 


11 36 47.97 


+54 17 37.25 


SAB(rs)bc 


21.7 


10.38 


[5] 


F658N (A) 


NGC 4041 


12 02 12.17 


+62 08 14.20 


SA(rs)bc? 


22.7 


10.55 


[2] 


F658N (A) 


NGC 6217 


16 32 39.22 


+ 78 11 53.60 


(R)SB(rs)bc 


18.3 


10.46 


[2] 


F658N(A) 



Note. — Coordinates (J2000) correspond to the UVIS aperture positio ns. Morphologies are tatc en fron i[de_ Vaucon leurs_et^jal.lj[T 9911 , 
in the NED homogenized notat ion), while d istances are drawn from: [ll lWilliclc et al.l | [T997|) : f2l lTullvl l fl988f) : [31 ITullv et al.l l l20od) : 
[4] iKe nnicutt et al. (2011); [51 ISpringob et a l. (2009). Stellar masses are derived through the 'best' Ks photometry from 2MASS 
ilSkrutskie et al.li2006i') . The final column ('Other HST data') lists non-proprietary imaging datasets from ACS (A), WFPC2 (W2), 
and NICMOS (N) in the MAST archive using filters other than B V7-equivalents. NGC 6217 was scheduled and observed with 
UVIS/F336W, but the observation failed. 



2.1. NGC 4041: Observations and Data 
Reduction 

Our chosen proof-of-concept object was NGC 4041, 
a face-on SAbc galaxy near the upper end of our 
distance hmit at 25 Mpc. This choice aUows for 
a demonstration of the methodology of source 
selection and processing for approximately half 
the SHUCS galaxies, those at distances beyond 
sa 15 Mpc. The processing and analysis of nearby 
systems will be outlined in a future paper. In 
addition, it offers an interesting study of envi- 
ronmental effects on star cluster formation and 
evolution in a structurally segmented system. 

We discuss the structure and environment of 
this galaxy in Section 2] and show a compos- 
ite HST image in Figure [H The HST observa- 
tions of NGC 4041 are shown in Table El The 
archival ^450, Vgoe, .^814, and Ha (F658N) data 
were obtained by two separate programs (9042 
and 9788). All data were retrieved from the 
Mikulski Archive for Space Telescopes (MAST) 
and each dataset was combined, corrected for 
geometric distortion, and drizz led to the native 
pixel scale using multi drizzle ( Koekemoeill2005t 
Fruchter fc Sosevll2009l) . For each WFPC2 image, 
two undithered exposures were obtained. These 
were drizzled separately for each CCD and the fi- 
nal pixel scale was 0"05 (PC) and O'.'IO (WF2, 3, 
4). The ACS F658N data consisted of a pair of 



images and these were drizzled to a pixel scale 
of 0'.'05. Only a single F814W ACS image was 
taken; cosmic rays were removed using the LA- 
COSMIC task ( van DokkumI l200l"l ) and the resul- 
tant image drizzled. For the WFC3/UVIS data, 
three dithered images were obtained and these 
were drizzled to a pixel scale of O'.'Ol. 

We have not corrected the images or photome- 
try for charge transfer efficiency (CTE) degrada- 
tion. The ACS images were taken two years af- 
ter the instrument was installed and CTE losses 



at th is epoch are negligible (jUbeda fc Anderson 
2012I ). CTE losses are apparent in the individ- 
ual WFC3/UVIS images but the background is 
above the critical level of 5 electrons where CTE 



losses are severe (see lNoeske et al.ll2012l ). Given 
that there was no correction available when we 
processed the WFC3 data, and the CTE losses 
are a small effect, we have not taken CTE losses 
into account. For WFPC2 data, no CTE trails 
are apparent in the images because of the fairly 
high backgrounds and we have therefore not ap- 
plied any photometric corrections to the WFPC2 
data. 

3. Selection of Star Cluster Candidates 
and Data Products 

Star clusters are only marginally resolved at the 
distance to NGC 4041, as the typical diameter of 
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Table 2 

Archival and new observations of NGC 4041 



Instrument 


Date 


Filter 


Exposure 
(b) 


Program ID 


P.I. Name 


WFC3/UVIS 


2011 Jan 30 


F336W 


1800 


12229 


Smith 


WFPC2 


2001 Jul 04 


F450W 


320 


9042 


Smartt 


WFPC2 


2001 Jul 04 


F606W 


320 


9042 


Smartt 


WFPC2 


2001 Jul 04 


F814W 


320 


9042 


Smartt 


ACS/WFC 


2004 May 29 


F814W 


120 


9788 


Ho 


ACS/WFC 


2004 May 29 


F658N 


700 


9788 


Ho 



7 pc subtends an angle of O'.'l on the sky. This 
translates to one or two pixels on various HST 
cameras. A considerable drawback is hence the 
potential inclusion of stellar associations, which 
are not discriminated by automatic source de- 
tection algorithms (see Silva- Villa fc Lars"enll20ld 



Bastian et al. I l201ll ). Toward that end we employ 



the concentration index (CI), defined as the dif- 
ference in brightness between two apertures: one 
comparable to the size of the point-spread func- 
tion (PSF), an d another representing the typical 
star cluster (see Whitmore et al.l[2007 ) . This tech- 
nique helps to place clusters between stars, which 
have very small CI, and associations, at large CI. 
We will develop the use of this method in and Sec- 
tion [Ol 



With the above in mind, we developed a selec- 
tion and photometry pipeline, which is run entirely 
in IRAELO and consists of the following steps: 

1. Source selection in [/336. We run DADf ind 

to select sources brighter than 7 a above 
background, measured on various parts of 
the image. We do not restrict 'roundness' 
(axial ratio) or 'sharpness' (size compared 
to the stellar full width at half-maximum), 
in order to include elliptical and marginally 
resolved clusters (the majority at this dis- 
tance, as noted above). 

2. Coordinate transformation. Our dataset 
consists of images taken for different projects 



t IRAF is distributed by the National Optical Astron- 
omy Observatories, which are operated by the Association 
of Universities for Research in Astronomy, Inc., under co- 
operative agreement with the National Science Foundation. 



at various times and with different cameras, 
which often leads to overlapping, but not 
identical pointings. Our U^sq observations 
were designed for maximum overlap with the 
archival imaging, and hence serve as the ref- 
erence frame for the World Coordinate Sys- 
tem (WCS). The coordinate lists from StepH] 
are converted from the frame to those of 
each instrument/pointing used in the study. 
This step uses the tmatch algorithm with 
typically some 20 reference sources (stars or 
compact clusters) that span the entire image 
as much as possible. The process is refined 
until the root-mean-squared errors for x and 
y positions do not exceed 0.1 px, a scale 
tested to eliminate source confusion. When 
registering WFPC2 data, the four CCDs are 
mapped individually, and the number of ref- 
erence stars is usually ten or lower. 

Multi-band photometry. Photometric 
apertures of radius O'.'l 2 are placed at the 
coordinates defined in Step [21 with the back- 
ground measured locally in annuli of 0'.'2 to 
0'.'3, depending on the pixel scale of each 
camera and the corresponding PSF. These 
values were derived after testing the typi- 
cal growth curves of mo deled stellar PSF s 
generated with Tinytim (jKrist et al.ll2011l ). 



coupled with the STScI Focus Model Util- 
ity^ to account for breathing at the time 
of observation. We note that these val- 
ues will be revised for each galaxy depend- 
ing on its distance. We do not allow the 
DAOphot task to re-center sources as the 



t http: / / www.stsci.cdu/hst / observatory/focus/FocusModel 
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differing pixel scale between cameras intro- 
duces non-negligible positional uncertain- 
ties. Given the employed pixel scales of 
the WFC3/UVIS (0'.'040 px"!), ACS/WFC 
(O'.'OSO px-i), WFPC2/PC (0'.'045 px'^), 
WFPC2/WF (Cy.'lOO px-i), the radius of the 
applied photometric aperture corresponds to 
image sizes of 3.0, 2.4, 2.7, and 1.2 px respec- 
tively. In the test case of NGC 4041, these 
translate to a physical radius of 13.4 pc, 
cf. the typical effective cluster radius of 
~ 3.5 pc. The photometric error is com- 
puted according to the DADphot recipe that 
takes into account the gain and readout noise 
of each detector. For cases where imaging 
from more than one HST camera is available 
in a given filter, we photometer all images 
and choose the measurement with the lowest 
associated error. We calibrate the photome- 
try in the Vegamag system. 

The photometric process is completed by the 
measurement of a CI for each source as the 
difference in C/336 brightness between aper- 
tures of one and three pixels. Such a tech- 
nique has been demonstrated in the past to 
produce viable samples of star clusters at 
comparable distances (e. g. The Antenna, 



Whitmore et al.l 120071 ). More information 



can be found in Section 

Photometric corrections. At this stage 
we correct for the size of the aperture. 
This has to be tailored to each galaxy in- 
dividually, as the correction is dependent 
on distance and whether the sources are 
resolved or not. In addition, we do not 
always have enough bright, isolated clus- 
ters to derive growth curves. The correc- 
tions thus combine our own empirical de- 
ductions with the encircled energy curves 
given for each instrument. For the test case 
of NGC 4041, isolated sources were very 
rare, so corrections were derived through 
the growth curves of model PSFs, gener- 
ated for each detector-filter combination 
with Tinytim. Finally, we add to our 
photometric errors a factor of 0.05 mag in 
quadrature, to account for uncertainties 
in t he photometric z eropoints (for details 
see lAdamo et al.l l2010 a| ). and correct fo r 



3.1. Completeness estimation 

At the distance to NGC 4041, the bulk of its 
cluster population is nearly point-like. There- 
fore, we follow the traditional method of artifi- 
cial star counts to estimate completeness. In fu- 
ture works focussed on nearby systems, however, 
we will use artificial clusters instead, to overcome 
the uncertainties introduced by stochastic vari- 
ations of the colors of low-mass clusters (intro- 
du ced by under-sampled stellar mass functions; 



see Fquesneau fc LanconllioiOl: Popescu fc HansonI 



2OIOI : ISilva- Villa fc Larsenll201ll ). 



^ 0.6 



r, 0.4 



E 
o 
o 




foreground extinction (jSchlegel et al.lll998[ ). 



22 24 26 

^336 [mag] 



Fig. 2. — Completeness tests tracing the fraction 
of artificial stars recovered at a range of bright- 
ness. The dotted vertical lines indicate 50% and 
90% completeness fractions for the inner disk and 
outer galaxy, performed separately to account for 
the different background levels. The top axis dis- 
plays the range of absolute magnitudes covered. 
We estimate 50% and 90% completeness levels at 
(25.5, 23.6) mag and (26.3, 24.6) mag for the inner 
disk and outer galaxy respectively. 



We created sets of artificial stars in a 15 x 15 ob- 
j ect g rid, using the MKSynth algorithm of LarsenI 
(|1999|). We used TinyTim models of the WFC3 
PSF as described in Step ^ of the pipeline (Sec- 
tion [S]). The generated objects were assigned mag- 
nitudes in the range [20, 30] mag, in steps of one 
magnitude. Since co mpleteness varies with the lo- 
cal background (e. g. IScheepmaker et al. 2007), as 
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well as crowding and confusion, we use two fields 
for this process: one covering most of the inner 
disk, and another, equally-sized field that covers 
the outer spiral structure. Results are presented in 
Fig. [21 where we show the 50% and 90% recovery 
fractions to occur at (25.5, 23.6) mag and (26.3, 
24.6) mag for the inner disk and outer galaxy re- 
spectively. 
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Fig. 3. — Brightness-error plots across the UBVI 
filter-set. The separation into inner and outer 
galaxy (blue dots, red triangles) accentuates the 
issues faced with photometric errors in the inner 
disk, an effect of crowding. This is taken into ac- 
count in all analysis. Using this plot we tentatively 
establish Veoe as the limiting band for this study, 
as it shows the most irregular brightness-error dis- 
tribution. Regular statistical measures, such as 
the standard deviation of the error distribution, 
provide no distinction between the BVI filters. 



Another photometric property that will impact 
our star cluster analysis is the limiting filter. We 
estimate this through the brightness-error plot of 
Figure [31 Regular statistical measures, such as 
the standard deviation of the error distribution, 
provide no distinction between the BVI filters. We 
therefore look for the most irregular distribution 
and adopt Vgoe as the limiting band. 



3.2. Concentration Index 

In SteplHof the photometry pipeline (Section[3]) 
we note the measurement of the CI, as the dif- 
ference in magnitude between apertures of ra- 
dius 1 and 3 px. We employ this metric in or- 
der to distinguish between star clusters and stel- 
lar associations, large , short - lived structures (see 



Portegies Zwart et al. 20ld Bastian et al. 2011 



20121) that can contaminate samples and affect the 



analysis of star cluster populations. 

In order to employ a CI cut in our source fil- 
tering process, we test the theoretical expectation 
from modeled clusters of various ages and sizes. 
Fir st we use MKSynth ( see Section IXTj) to convolve 
an lElson et al.l ^9871) surface brightness profile 
with a PSF derived from point sources in our im- 
ages. Since we could not identify suitable stars 
on the F336W image of NGC 4041, we employed 
SHUCS images of NGC 891 and NGC 2146. In ad- 
dition, while various images are taken at different 
focus settings, the PSF does not change to a level 
that will affect our analysis - the effect is smaller 
than our typical photometric error of 0.1 mag. 
We assign a range of sizes to these model clusters, 
incorporate them on the F336W image and mea- 
sure their CI in the same way as for the detected 
sources. The top panel of Figure [H shows a plot of 
CI versus assigned radius. The first datapoint, at 
-Rcff — pc, denotes individual stars, which reg- 
ister values of w 1.0 mag. We therefore adopt a 
lower limit for cluster candidates at CI = 1.1 mag. 
We adopt an upper limit from the literature, as 
the CI of the largest observed clusters of approx 
imatel y 10 pc (see si z e dis tr ibutions of iLarsen 
2004bt iBarmbv et all 120061: [Scheepmaker et al ' 



2007l : lMavva et al.ll2008llBastian et al.ll2012l) . The 

corresponding CI is 1.8 mag on Figure [H 

To confirm the asymptotic shape of this relation 
we employ clusters generated as aggregates of indi- 
vidually synthesized stars, fo l lowing the method- 
ology of ISilva- Villa fc Larsen ( 2010t ). Each star is 
represented by the empirical PSF described in Sec- 
tion l3.1l and a King profile of index c = -tidal _ 3Q 
while the age, mass and effective (half-light) radius 
assume the following ranges: 

log(A//M0) = [3.0,3.5,4.0,4.5,5.0] 

log(T/yr) = [6.6,7.0,8.0,9.0] 
i?cff= [1,2,3,4,5,7,9,12,15,20] pc 



8 



The range of CI assigned to clusters is shown as 
a hashed area in the top panel of Figure HI while 
in the center and bottom panels, the CI area is 
bracketed by dotted lines to aid illustration. The 
CIs of stochastically modeled clusters are slightly 
higher (by « 0.1 mag), owing to their King sur- 
face brightness profiles, chosen to generate more 
realistic older clusters. 

The measured CI values of actual sources are 
plotted against C/336 brightness in Figure [5j Ap- 
plying the above CI cut and the 90% completeness 
limits of Section IO] we restrict cluster candidates 
to a very specific part of photometric parameter 
space. 

3.3. SED-fitting with Yggdrasil models: 
derivation of age, mass and extinction 

The age, mass, and extinction of all candi- 
date star clusters are derived through fits to their 
spectral ene rgy distribution (SEP) with Yggdrasil 
models (see IZackrisson et al. 2011 . for a detailed 
description). Such a model set is ideal for our 
work, as it combines single stellar population 
(SSP) synthesis models with nebular spectra, both 
emission lines and continuum, and is tuned to 
the redshift of each galaxy. Analyses using Yg- 
gdrasil models have in the past been successful 
in describing the cluster p opulations of intensel y 
star forming galaxies (e.g. Adamo et al.l bOlObh . 
For the SHUCS sample, we use Starburst99 



SSPs (jLeitherer et al.L l999: Vazquez fc Leitherer 
2005 ). compiled using a iKroupal ( 2001 ) stellar 
initial mass function throughout the mass range 

0. 1 — 100 Mq; Padova stellar evolutionary tracks; 
and three different metallicities: 0.4, 1, and 2.5 
times solar. We note no spatial trends in terms of 
metallicity, as shown in FigurelHl While this result 
might seem to oppose the well-established metal- 
licity gradients in local spirals (e. g. in the Local 
Group; |Cioni 2009) , the photometric derivation of 
metallicity is subject to large uncertainties. No 
firm conclusions can be drawn about the metal- 
licity distribution. To reproduce the integrated 
fluxes of the earlier stages of cluster evolution, 

1. e. luminous H il regions surrounding the recently 
formed stars, the models include a self-consistent 
treatment of the nebular emission p erformed with 



the C loudy photoionization code (jFerland et al 
19981 ). The metallicities of the gas and the stars 
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Fig. 4. — Concentration Index tests on artificial 
clusters. The top panel s hows a CI v e rsus r adius 
plot for clusters created as lElson et al. ( 1987 ) light 
profiles of index 7 — 3.0, convolved with an empir- 
ical PSF. The center and bottom panels show ar- 
tificial clusters created as aggregates of individual 
stars, arranged according to a King profile (more 
appropriate for older clusters). Cluster ages are 
log(r/yr) = 7, 8, while different symbols represent 
a range of masses. As they age, clusters display 
less regular behavior, as evidenced by the error 
bars in the log(r/yr) = 8 model set. A consistent 
asymptotic behavior is noted through both meth- 
ods. Rctf > 10 pc are rarely observed, so we set the 
CI at this radius as the upper limit. An Res = 
represents individual stars, and the corresponding 
CI ~ 1.0 mag is adopted as the lower limit. 



are assumed to be the same. Using the gas proper- 
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Fig. 5. — Concentration Index versus source 
brightness on isotropic axes. The bottom panel 
shows all sources that register UBVI photometry 
as green points, while the open and closed gray 
circles apply a cut according to photometric error. 
The top panel shows the space chosen to represent 
high-quahty cluster candidates, the ones used for 
all further analysis, against a backdrop of the gray 
circles of the lower panel. Vertical lines denote the 
90% completeness limits for the inner and outer 
galaxy. 



tiesof typical local H ii regions (jKewlev &: Dopita 
2OO2I ) we assume a filling factor / = 0.01, a hydro- 
gen density of n(H) = 100 cm~^, and a covering 
factor, i. e. fraction of ionizing photons absorbed 
by the gas, of c = 0.5. Different values of c, 
n(H), and / only impact the SED fits for the very 
youngest clusters. We therefore never attempt to 
resolve the 1 — 5 Myr age block, and accept a 
single, 'young' age for such sources. 

These models, together with the photometric 
catalogs, comprise the i nput for a least-y^ fit ting 
algorithm, as detailed in lAdamo et al.l (j2010a[) . In 
brief, the program considers sourc es detected in 



four o r more filters and applies the iCardelli et al 
attenuation law at each age step to con- 
strain the visual extinction in the line of sight. 
The final age, mass, and extinction are assigned 
to a cluster given the set of parameters that min- 
imizes the x^- Uncertainties are carried through 
the fitting process. The program estimates the 
best reduced xt best (x^ divided by the degrees of 



freedom), and saves, as the range of acceptable val- 
ues, all the solutions with < xtbest + The 
maximum and minimum age, mass, and extinc- 
tion contained in this pool of likely solutions are 
assigned as errors. This is a typical method of es- 
timating errors when the statistic i s empl oyed, 
e.g. the photometric fits of Bik et al] ( 20031) and 
the sp e ctroscopic fits of Kqnstantopoulos et al.l 



(200' 



Lampton et al.l (|1976[ ) provide a full rea 



sonmg. 

Example SED fits are shown in Figure [T] First 
we provide an illustration of the best fitting Yg- 
gdrasil model, with source photometry marked as 
filled dots. The red squares mark the best-fit 
model fiux on the y axis, along with the central 
wavelength and bandpass of each employed HST 
filter. The plots of the bottom row plot in 
the age-extinction and age-mass spaces. Red con- 
tours mark the range of acceptable solutions to 
the SED fit, defined as xl < xl,best + ¥■ The top 
row of figure sets features two cluster candidates 
that meet the selection criteria (detailed in next 
section), while the ones in the bottom row do not. 

3.4. Cluster candidate selection and data 
products 

The primary data products of each study are 
a Photometry Table and an SED Table. These 
data products will be made available to the com- 
munity in their entirety when the survey has been 
completed. The Photometry Table collects posi- 
tional and photometric measurements and applies 
a photometric flag, fuBVi, according to the follow- 
ing scheme: 

1: Source detected and photometered in UBVI, 
11: ... plus UBVI errors less than 0.3 mag, 
111: . . . plus a CI in the range [1.1, 1.8] mag, 
1111: ...plus brighter than the 90% com- 

pleteness limit. 

All analysis is performed on clusters with 
fuBVi — 1111. The scheme is complemented by a 
second flag, fuBVi+, which also considers the pho- 
tometric quality in bands other than UBVI (in 
this case Ha), and is marked with a similar se- 
quence of the number 2 for clarity (2, 22, 222, 
2222). The effect of filtering the sample across 
the UBVI baseline on photometric completeness 
is negligible, and will be discussed fully in the 
context of the luminosity functions of Section [5.2l 
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Fig. 6. — Metallicity distribution (left) of all sources in the filtered catalogue {fuBVi — 1111, see Section [3^ . 
divided according to galactocentric radius. We find sub- and super-solar metallicities to be preferred in 
NGC 4041. However, the uncertainties in photometrically derived metallicity are considerable, so no firm 
conclusions can be drawn from this histogram. The metallicity distribution also does not show any spatial 
trends. This is shown on the right, where all clusters in the final sample are symbol coded according to 
derived metal content. 



This filtering process, as applies to NGC 4041, 
is demonstrated in Fig. [S] The top row shows 
two segments of our F336W image where we mark 
sources with fuBVi of 1 through 1111 as red 
dots and green, blue, and yellow circles. The ra- 
dius of 3 px matches the photometric aperture in 
17336. The bottom row shows a UBVI color di- 
agram, color-coded in the same way, to demon- 
strate the effect of each filter on the selection 
process, along with the corresponding brightness- 
error plot. From these images it is clear that the 
distance to this galaxy gives rise to a crowded field 
of clusters and associations. Coupled with the 
stochastic variations in the light profile (present 
at all distances) , this severely cuts down the num- 
ber of clusters for which reliable size measurements 
can be performed. At 22.7 Mpc one UVIS pixel 
corresponds to w 4 pc, while one WFPC2/WF 
pixel is 11 pc. PSF deconvolution codes (e.g. 
IShape; iLarsen Il999l) have been demonstrated to 
resolve features as small as 0.4 px on HST imag- 
ing, which corresponds to just under 2 pc. That 
way, the smallest effective radius that can be sam- 
pled is w 1 pc. It is doubtful that such a source 



would have a sufficiently high signal-to-noise ratio 
to be detected, however, as a test run of IShape 
on NGC 4041 fit only large clusters, typically 
ReS > 5 pc. This topic will be revisited in future 
SHUCS works that treat nearby systems. 

In addition to the photometry flags, the Pho- 
tometry Table also marks each source with an in- 
teger between 1 and 4 to denote which WFPC2 
CCD covers the source. A value of '9' marks a de- 
tector other than WFPC2, while '0' indicates that 
a source is covered by WFPC2, but it is too close 
to the edge of the image. This flag is then repeated 
for each filter in the analysis, standardized as U, 
B, V, Ha, I, J (e.g. NICMOS-FllOW), H (e.g. 
NICMOS-FllOW). The shorthand filter represen- 
tations are defined in the table header. Positional 
information is also included, in the form of coor- 
dinates measured on the F336W image. 

The results of the SED fits are collected in an 
SED Table, featuring a minimum, best-fit, and 
maximum value for age, mass, and extinction; the 
number of data points (bands) employed by the 
fit; and the reduced statistic of each fit. 
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Fig. 7. — SEP fits of four cluster candidates in NGC 4041 with Yggdrasil models, following the process of 
Adamo et alj (|2010al) . The filled blue dots represent the observed integrated fluxes, with associated error 
bars. The corresponding integrated model fluxes are labelled as red squares and the horizontal lines indicate 
the bandpass of the given HST filter. The quality of the fit is shown in the lower panels of each set, where 
the is plotted across the parameter spaces of age-extinction and age-mass. The red contour line indicates 
the range of acceptable fit solutions. Each plot is labelled with an identifier, and the derived age and mass 
of the cluster candidate. The top row shows good fits, while the cluster candidates of the bottom row are 
filtered out of the final catalog. Cluster 4427 is found in a crowded region and its light is contaminated by 
neighboring sources. As such it does not pass the CI test employed by the filtering process and receives 
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Fig. 8. — The cluster selection process, demonstrated for the population of NGC 4041. The top panels show 
close-ups of the F336W frame (inner, outer galaxy on the left and right). Each image measures at approxi- 
mately 1.3 X 0.9 kpc. We mark fuBVi values of 1, 11, 111, and 1111 as red dots and green, blue, and yellow 
circles of radius 3 px, coinciding with the photometric aperture. The bottom panels follow this select ion on 
a UBVI color diagram, with associated single stellar population model tracks (jZackrisson et al. I l201l[ ). The 
metallicities of the dashed (green), solid (black) and dashed-dotted (magenta) lines are 0.4, 1.0, 2AxZq. On 
the right we show a J/sae magnitude-error plot. None of the criteria appear to introduce color-space biases. 



The Structure 
NGC 4041 



and Environment of 



(jl988[ ). Given this distance modulus, we derive 
an absolute magnitude of 5 = —19.7 mag, accord- 
ing to i ts Sloan Digital Sky Survey Data Release 6 



mg to I ts bloan Uigitai bRy burvey Data Kelease b 
We begin the analysis of the SHUCS dataset iSDSS: lYOTiret al.l2000HAdelman-McCarthv et al 



with the SA(rs)bc galaxy NGC 4041 ()de Vaucouleurs et al. 



19911). We adopt a distance of 22.7 Mpc and a cor- 



20081 ) model or Petrosian radius brightness. 



responding distance modulus of 31.78 from iTuUv 



In Figure[T]we showed an HST color-composite 
image of the galaxy in UBVI and Ha. The 
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galaxy appears to have a partitioned structure, 
with tightly wound inner spiral arms giving way to 
a looser, outer network. The boundary is defined 
by a smooth, yellow bulge-like feature, suggestive 
of an old underlying stellar population. The bulge 
is overlaid with bright, blue star clusters and dust 
lanes, and it fades quite rapidly toward the out- 
skirts of the galaxy, where blue star clusters can 
be found in large numbers. Finally, we observe 
the nuclear source to be split into two components 
with comparable fluxes, both offset from the center 
of the outermost isophote by O'.'S and 1". The fol- 
lowing sections treat the galaxy and its surround- 
ings in more detail. 

4.1. A complex isophotal structure 

Multi-component spirals have been observed in 
the past, first th e southern supergiant NGC 6902 
([Gallagher Il979l) . and then in the context of low 
surfa ce brightness galaxies (Bosma fc Freeman 



Galactocentric radius rkpcl 
2 4 6 B 
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19931 ) ■ 

such structures 



Matthews fc Gallagheij (|l997t ) observed 
in a sub-sample of 'extreme' 



late-type galaxies. They also find offset cen 



tral p eaks in such galaxies, as do lOdewahn et al 
(Il996) in another catalog of late- types, and 



Karachentsev et al. (Il993 ) in Magellanic spirals. 
Matthews fc Gallagheil discuss the possibility of 
offset central peaks representing the true dynam- 
ical centers of their hosts, i. e. centers of grav- 
ity offset from the ge ometric or isophotal center. 
Marconi et al. (l2003l) touch on this scenario by 
suggesting that the dynamical center of NGC 4041 
might be decoupled from its geometric center. 

The composite structure of NGC 4041 is investi- 
gated further through plotting the surface bright- 
ness of the galaxy in azimuthally averaged bins, 
normalized to the central peak (despite the small 
central offset, to simplify illustration). This is 
shown in Figure [9l where the central exponen- 
tial disk measures at a radius of « 22" (2.4 kpc) 
and we derive a radius of ~ 80" (9 kpc) for the 
outer galaxy from the optical isophotes. As noted 
above, the spiral structure is not discontinuous, 
with some spiral arms stemming from the floc- 
culent inner structure and developing into loose 
outer arms. Across the disk the inter-arm dis- 
tance increases from ^ 0.7 kpc to « 1.1 kpc, as 
measured on the F336W image, and the arm thick- 
ness is « 0.5 kpc. The inner disk follows an ex- 
ponential profile, as is expected for the disks of 
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Fig. 9. — Peak- normalised, azimuthally-averaged 
count rate of NGC 4041 in bins of 3'.'5 (400 pc) 
about the central peak. While the inner disk is 
adequately described by exponential flux profiles, 
the outer structure follows a profile more suited 
to early- type gal axies - a Sersic index of 4.0 is 
equivalent to the de Vaucouleurd ( 19481 ) r^/^ law, 
while an index of 1 equals an exponential profile. 



spiral gal axies, while t he outer disk is better de- 
scribed bv lSersic ( 19681 ) profiles of high index, typ - 



ical of early-type galaxies (e. g. Cote et al. 2006[ ) 



The F336W profile shows a periodic fluctuation 
While one might expect the popu lation of youuj 
star clusters to contribute to this, Larsen ( 2004 



show the contrary: star clusters are never found to 
contribute more than 10% of the overall U band 
luminosity of a spiral galaxy (and up to 20% in 
starbursts). 

4.2. Neighbors and Potential Companions 

The complex structure of NGC 4041 could 
be the result of a past dynamical event, so we 
search the NASA Extragalactic Database (NED) 
for neighbors at accordant redshifts. We limit 
the search to galaxies in a cone of 1 degree 
about NGC 4041 and in a velocity space of 
Av < 400 kms^^. The search yields five neigh- 
bors, so the galaxy is by no means isolated. Details 
are given in Table [H and a finding chart shown in 
Figure [lOl 

The first five galaxies in t his list c o mpris e group 
LGG 266 in the listing of ICarcial (|l993l) . The 

sixth, MGC-Ml-15-013, also fits the LGG cata- 
log limiting criteria - a velocity difference of less 
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Fig. 10. — The loose grouping of galaxies around NGC 4041. One is a distorted spiral, three appear irregular 
(albeit the high inclination hinders this classification), and one is an early- type galaxy with a faint inner 
dust lane (evident in regular contrast scaling, but not here). We mark the velocity of each galaxy on the 
SDSS mosaic image on the left, and show high-contrast zoom- ins on the right (not to scale). The four 
smaller galaxies (top two rows) also exhibit considerable FUV fluxes on GALEX images, indicative of recent 
or ongoing star formation. Their calculated SFRs are listed in Table [3] The color of the inner disk of 
NGC 4041 appears similar to that of the lenticular NGC 4036, while the outer disk shines bright in shorter 
wavelengths. 



than 600 kms~^ and a projected separation less 
than 0.52 Mpc - but was only discovered recently 
( Liske et al Since we find no morphological 

classification for this galaxy in the literature, we 
inspect its z band SDSS image to find a linear pro- 
file with no pronounced central peak. We therefore 
classify it as Irr. The close-up of Figure [TUl shows 
a clumpy profile, but it is biased by bright star 
forming regions at the extremes of the disk that 
are not detectable in the z band. The dominant 



galaxy of this group is NGC 4036, a lenticular. It 
features an inner equatorial dust lane and registers 
LINER emission (Veron-Cett v &: Veron.2006i) . 

4.3. Extended Structure 

Upon inspecting GALEX UV imaging of 
NGC 4041 we find a faint, three-pronged feature 
to the south of the UV disk. It extends between 
125" and 220", or - 14 - 24 kpc. This 'tidal 
fork' is shown in the far-ultraviolet plus near- 
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Table 3 
The NGC 4036 Galaxy Group 



ID 


RA 


DEC 


rUg 


VR 


SFR 


log(M.) 


Morphology 


References 




(h m s) 


(o / //) 


(mag) 


(kms-^) 


(Mo yr-i) 


(Mo) 






NGC 4036 


12 01 26.7 


+61 53 45 


11.97 


1385 


0.02 


10.42 


SO 


[1] [1] [2] 


NGC 4041 


12 02 12.2 


+62 08 14 


12.19 


1234 


0.63 


10.58 


SA(rs)bc? 


[1] [1] [3] 


IC 0758 


12 04 11.9 


+62 30 19 


14.25 


1275 


0.19 


9.26 


SB(rs)cd? 


[4] [4] [3] 


UGC 7009 


12 01 44.1 


+62 19 33 


14.33 


1119 


0.12 


8.78 


Im 


[4] [4] [4] 


UGC 7019 


12 02 29.4 


+62 25 02 


15.05 


1518 


0.05 




Im 


[5] [6] [7] 


MGC+11-15-013 


12 02 43.3 


+62 29 52 


16.07 


1059 


0.03 




Irr 


[4] [4] [4] 



Note. — Morphologies are taken from lde Vaucouleurs et aP l)l99ll l. apart fr om MGC+11-15-013, which we classify in 
SectionSm Stellar masses are derived though 'best' 2MASS Ks magnitudes (jSkrutskie et al.ll2006l l. except NGC 4041, 
for which we use an aperture of 72'.' to match the measured aperture of the outer disk. SFRs are derived through 
GALEX FUV imaging (see Section I4.4[l . SDSS photometry is model photometry in the g band, where applicable. 
Positions, photo metry, and radial velocities are dr awn from the following sour ces, quoted in triplets in co l umn 7 : [1] SDSS 
Data Release 6. lAdelman-McCarthv erall (l2008h: [2] ICappellari etHI (|201ll ): [3] Ide Vaucouleurs et all (Il99ll): [4 ] SDS S 
D ata Release 2, lAbazaiian et all (|2004h : [5] ICotton et al.l (|l999l ): [6] VATT B band photometrv. iTavlor et al.l (l2005h : 
[7] ISpringob et al.l ( 20051 ): The first five galaxies in this list make up group LGG 266 in the listing of Garcial 1 19931 ). and 
the sixth also obeys the inclusion criteria. Given the individual distance moduli, we derive absolute g band magnitudes 
of -20.0, -19.6, -17.9, -17.3, -16.7, and -15.0 mag. 



ultraviolet (FUV+NUV) image of Figure [TTl We 
stack the two GALEX channels to boost the flux, 
and smooth by a three-pixel Gaussian kernel to 
remove noise. We present this image at high con- 
trast to accentuate low surface brightness features. 
The contours mark optical brightness, derived by 
stacking the g and r SDSS images, and smoothing 
by a ten-pixel Gaussian kernel. 

The feature is detectable in the UV image at 
the 3-5 a level, but not quite as bright in the 
stacked and smoothed optical image (< 2 cr de- 
tection). Assuming this is a tidal feature physi- 
cally associated with NGC 4041, we can attribute 
it to a past interaction. The low 5 + r flux denotes 
low or no Ha emission, and therefore no current 
star formation in the 'fork'. The UV flux should 
then mostly originate from activity on the order of 
100 Myr ago, rather than current star formation 
and 0-stars. 

4.4. Star Formation Rate, Stellar and Gas 
Mass 

Since there are no published values for the 
star formatio n rate of NGC 4041, we apply the 
Kewlev et al.l (|2002( ) methodology to obtain esti- 



mates from the IRAS infrared fluxes a t 60 //m and 



lOOum. In brief, we flrst employ the iHelou et al 



(|1988[ ) prescription to obtain the far-infrared flux, 
FpiR, from the IRAS fluxes, which we then con- 
vert to a luminosity through Lfir = 4 tt x i^FiR, 
where D is the dista nce to the galaxy . We then 
use th e lKewlev et al.l adaptation of the lKennicutt 
(jl998l ) law to derive a FIR star formation rate 
(SFR) of 4.10 Mq yr-i for NGC 4041. This 
compares well with galaxies of similar morpho- 
l ogical type and brigh tness in the SINGS survey 
( Kennicutt et al.ir2003[) . The flux ratio between 60 
and 100 ^m of 0.44 is also in accord with SINGS. 

We also obtained an FUV SFR for NGC 4041 
through archival GALEX imaging, which we ex- 
pect to be lower, due to the absorption of UV 
photons by dust. To estimate the flux from 
the target we flt it with the ELLIPSE function 
in IRAF and added up the counts correspond- 
ing to the two segments of the galaxy disk (as- 
suming the above radii of 80", 22" and no el- 
lipticity), while accounting for a flat background 
level and a foreground extinction of E{B — V) — 
0.15 mag (following Ap uv = 7.9 x E{B - V) from 



Gil de Paz et al.l 120071 ). Note that the low reso- 
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is common in post-i nteraction systems, where gas 



Fig. 11.— FUV+NUV image from GALEX, 
smoothed with a three-pixel Gaussian kernel, and 
presented at high contrast to emphasize low sur- 
face brightness features. We reveal the presence 
of material beyond even the extended UV disk of 
NGC 4041, in the form of a three-pronged struc- 
ture. This 'tidal fork' might be part of a greater 
network of debris, but it is the only feature de- 
tected at a level of 3-5 a above the background. 
The contours show the isophotal structure of a 
g + r image from SDSS at the 1, 2, 3, 5, 10 a level. 
Only the brightest, central prong is detectable in 
SDSS at the 1-2 a level. The absence of Ha emis- 
sion (transmitted in the r band filter) suggests an 
age on the order of 100 Myr. Concentric circles 
mark the outer disk of NGC 4041, a faint, ex- 
tended UV disk, and the full radial extent of the 
fork. The angular distances correspond to 9, 14, 
and 24 kpc at the distance of 22.7 Mpc. The bright 
feature in the southwest is a foreground F-type 
star. 



lution of IRAS images does not allow for a sim- 
ilar treatment. We then convert fluxes to star 
formation rates of . 18, 0. 65 M0 yr~^ through 
the FUV iKennicuttI (|l998[ ) relation. Consider- 
ing the area of each segment, we derive SFR 



densities of [0.45, 0.08] M© yr ^ arcmin ^, or 
[1.03, 0.20] X 10"^ Mq yr~i kpc^^ for the inner 
disk and outer galaxy respectively. We apply the 
above methodology to also derive the FUV SFRs 
of the group around NGC 4041, and list results in 
Table |3l The observation of a markedly higher 
SFR density in the central regions of a galaxy 



20031: 



is tunneled inwards (iFqrster Schreiber et al 

Konstantopoulos et al .I l2009h until it is the nucleus 
alone that is experiencin g a starburst, or un til it 
ignites an active nucleus ( Ellison et al.ll201lh . 



Stellar masses, M ^, were derived by applying 



the iBell et al.l (|2003[ ) prescription with a Ks solar 
brightness of 3.32 mag to catalogu ed 2MASS Ks 



photometry ( Skrutskie et al. 2006|). The faintest 



two members of the group are below the detection 
limit, and the obtained values for the four brighter 
galaxies are listed in Table [3l While we generally 
use the 2MASS 'standard' radius, in the case of 
NGC 4041 we use a radius of 72" in order to en- 
velop the entirety of the outer galaxy. This gives 
rise to a larger for NGC 4041 than NGC 4036, 
despite the brighter tabulated g band magnitude 
of the latter. Additionally, we derive a specific 
SFR (sSFR) of NGC 4041, through division with 
the FIR SFR, of 10.9 x 10"" yr^\ consistent 
with the morphological type of NGC 4041 - cf. 
its near morphological counterparts HCG 7C and 
HCG 59A, that also r egister values ~ 10" yr"^ 
(|Konstantopoulos et al . 2010, 20il)- 

Finally, we look for information on t he gaseous 

compo nent of NGC 4041. The survey of lCouto da Silva fc de Freitas P 

(Il989l) offers a value of log(MHi/M0) = 9.76. 

We also use the C O luminosity published by 

Elfhag et al.l ( 19961 ) to estimate the mass of 

molecular hydrogen in the system. We use the 

classical Milky Way conversio n factor, aco = 



19871) . to convert the 
^ pc^ to a gas mass 
The combined gas 
Mgas = 1-3 Mhi + Mmoi (the factor of 



4.6 Mq ([Solomon et al 
log(Lmoi) = 8.48 K kms" 
of log(AW/MQ) = 9.14 

mass, 

0.3 M hi accounts for helium, e.g. iHunter et al 
Il982l ). is therefore log(Mgas/M0) = 9.95. Divided 
by the SFR, we get a gas depletion timescale of 
- 2.2 Gyr. 

4.5. A Possible Dynamical Event in the 
Recent Past 

Many of the traits examined in this section 
are consistent with a past dynamical event in 
NGC 4041: a segmented brightness profile; a dou- 
ble central peak, offset from the isophotal center; 
a UV-bright tidal feature; and a markedly higher 
central SFR. Of particular interest is the offset 
central peak, which potentially represents the dy- 
namically decoupled core of NGC 4041 (following 
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Marconi et al. 


2003). Such features herald past 


mergers (e. g. 


Emsellem et al. 2011). not interac- 



tions. 

The group environment is thought to be 
conducive to dynamical events and accelerated 
galaxy evolution 



Konstantopoulos et al. 20101 " 2012D . Furthermore, 
while a dense environment can enhance star for- 



.. , . , .e . g. iVerdes-Montenegro et al 

20011: IJohnson et all 120071: iGallagher et al.l bOloT 



mation during interactions (jMartig fc Bournaud 
20081 ). the effects of minor events are often not 
very pronounce d, leaving but the faintest tr aces 
of their passing ( Konstantopoulos et al. 2010[ ). At 
the same time, compact galaxy groups have been 
nomi nated as the sites of lenticular galaxy forma- 
tion (jWilman et al.ll2009l) . as minor mergers build 
up mass and gradually exhaust the gas reservoir of 
the group, while retaining the structure of individ- 
ual disks. The high-index Sersic outer brightness 
profile for NGC 4041 (Figure [3]) is more akin to a 
lenticular than a spiral. 

The overall symmetry of the disk seems to rule 
out a major merger. We therefore propose an ac- 
cretion event, perhaps that of a gas-rich dwarf, as 
the origin of the present star formation activity 
and the two-component disk of NGC 4041. 

5. The Bimodal Star Cluster Population 
of NGC 4041 

In the Introduction of Section [T] we argued that 
the understanding of a galaxy and its surround- 
ings can be enhanced by a study of the star clus- 
ter population. We now present such a study for 
NGC 4041. We split the sources according to 
galactocentric radius (rgc), in order to contrast 
the sub-populations of the inner and outer regions, 
and derive analysis only from high-confidence de- 
tections, i. e. sources flagged in Section 13.41 as 

fuBvi - nil. 

5.1. Star Cluster Colours 

The colors of star clusters are indicative of the 
age of the stellar population they represent, hence 
SSP models can be used to interpret the colors of 
cluster candidates across NGC 4041. The color- 
magnitude diagram (CMD) of Figure [12] plots the 
evolution of SSP models of various masses (from 
10^ to 10^ Mq), and shows dissimilar distributions 
for clusters in the inner and outer galaxy (blue 



circles/red triangles). While this hints at different 
mass functions, incompleteness will eliminate the 
low-luminosity end of the inner disk distribution. 
Still, such dramatic contrasts between subpopu- 
lations within a single system are not normally 
obse rved outside merging and interacting galaxies 
(e. g. lFedotov et al.l201l[ ). We remind that no spa- 
tial trends were discovered in the metallicity dis- 
tributions of clusters in the inner and outer galaxy 
(see Section 1331): and note that the reddening dis- 
tribution shows no change between inner disk and 
outer galaxy. 
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U336 - B450 [rnog] 

Fig. 12. — Color-magnitude diagrams of high- 
confidence cluster candidates with Yggdrasil mod- 
els of solar metallicity overplotted on isotropic 
axes. Model tracks represent masses of 10'^, 10**, 
10^, and 10^ M0 from bottom to top, each trac- 
ing ages between 6 Myr and ~ 10 Gyr. Num- 
bers in filled circles mark each age dex of model 
SSP evolution. There is a pronounced difference 
in the color distributions of clusters in the inner 
disk (blue dots) and outer galaxy (red triangles), 
likely a combination of differing detection limits 
and possibly also mass distributions. Such dispar- 
ities are often observed in interacting systems, but 
not elsewhere. The crosshairs in the bottom left 
represent the median error in each axis. 



We further investigate this disparity through 
the color-color plot of Figure [131 As in the CMD 
of Figure [T^l the model track spans the full evo- 
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lution of an SSP, from 1 Myr to about 10 Gyr. 
The SSPs account for the transmission of nebular 
continuum and emission hues in the F606W fiher. 
This extends the purely stellar SSP toward greener 
values and therefore cover sources that scatter to- 
ward the top left of the color distribution. We 
show color contours on the right panel to aid com- 
parison, where blue and red shades represent the 
inner/outer galaxy. To avoid over-representing 
outliers we omit the lowest two contour levels, 
hence plotting only two-dimensional bins that 
sample at least 15% of the full population. Sources 
in the inner disk appear to often diverge toward 
bluer values. We attribute this to crowding and 
aperture effects: on the vertical axis, the -8450 
flux is measured from WFPC2, at lower resolu- 
tion; and on the horizontal axis, we might expect 
some contamination by the diffuse, irregularly dis- 
tributed Ha emission in the inner disk. In addi- 
tion, the colors of low-mass clusters might pref- 
erentially diverge toward this region of colorspace 
(iFouesneau fc Lancon 2010 : Popescu fc Hanson 



I2OIOI : ISilva- Villa fc Larsenll201l[ ). A third possi 
bility is a variety of filling factors in the gas, which 
might lead to variations of up to 0.5 mag in the 
model tracks. Since we adopt a given filling factor 
and electron density, we might be underestimating 
the blue emission from these regions. 

The loci of the color contours (on the right), in- 
dicated by stars, are offset by w 0.5 mag, while the 
overall distributions extend each other: younger 
ages are under-represented in the inner disk, and 
the outer galaxy shows few clusters older than 
1 Gyr. This appearance is likely due to selection 
effects, namely the different detection limits in the 
inner disk and outer galaxy. On the young end, 
crowding in the inner disk might be responsible 
for this apparent lack of sources, especially ones 
with low mass. Furthermore, the higher SFR of 
the inner disk (Section 14. 4p is expected to corre- 
late with the number of young clusters. 

Such a disjointed color distribution might be 
expected in a mer ger, cf. the cluster population s 



of the Ant ennae ( Wh itmore & Sc hweizei 19951 ). 



NGC 3256 (|Trancho et al..,2007a. b'). and most no- 
tably NGC 7252, which exhibits an inner star 
form ing disk and an outer halo- type dis tribu- 



tion (jMiller et al.ll997l:ISchweizer fc Seit^ri998„ ). 
Even so, such an orderly inner- versus-outer galaxy 
distinction has not been noted in past analyses. 



5.2. Luminosity Function 

The Luminosity Function (LF) of star clus- 
ter populations is routinely observed to have a 
smooth, power-law shape, as noted in the Intro- 
duction. While this shape can be seen in the 
binned LF of Figure[T4](left), the cumulative func- 
tions on the right panel reveal a striking substruc- 
ture of: (i) an incomplete part with an evolving 
gradient; (ii) a smooth, power-law segment; and 
(iii) a truncation at the bright end. This indi- 
cates that substructure may have been habitually 
ignored in the literature by arranging clusters in 
bins of arbitrary width (see the cautionary note of 
Mai'z Apellaniz fc Ubedall2005l) . 

We measure the slope from the cumulative func- 
tion, in two ways. First we fit only the smooth 
part of each function, but maintain the fit consis- 
tent between the four bands. Then we fit over a 
large range, in order to emulate the effect of bin- 
ning. The results, listed in Tabled show shallower 
functions for the full range fit (consistent with the 
literature but slightly steeper than the canonical 
value of —2), while selective fits of the power- 
law segment return steeper values. Fitting only 
the power-law part of each function provides the 
added advantage of not biasing the measurement 
according to physical parameters such as cluster 
age or reddening. Since young clusters tend to 
be much brighter than older ones, the bright end 
of the LF will largely represent the young part of 
the population. Our tests indicate that fitting the 
entire range can lead to very large errors and pos- 
sible biases that are not present when fitting only 
the power-law segment. 

The C/336 returns markedly shallower fits than 
the rest of the bands, an eff ect noted in the liter- 
ature review of lGielesI (|201fl[ ). where redder bands 
were shown to have steeper distributions. We note 
that the extinction d oes not change the slope of 
the LF (Larsen 2002). The GielesI analysis also 
showed LF slopes to be a function of the mean lu- 
minosity. The mean brightness over which we fit 
(21 or 22 mag; see Table |4]) translates to a lumi- 
nosity of log(L/L0) ~ 5, corresponding to litera- 
ture slopes of —2.5. We attribute the consistently 
steeper slopes derived from the binned histograms 
to sampling uncertainties, especially given the rel- 
atively small cluster population being fit. This dis- 
crepancy exemplifies the loss of information that 
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Fig. 13. — Left: Usse — S450 vs Veoe — hu colors of star cluster candidates in NGC 4041 on isotropic 
axes. We overplot Yggdrasil models, following the conventions of Figure [T^l and also include sub- and 
super-solar tracks in this case. Clusters in the inner disk and outer galaxy are marked as blue circles and 
red triangles respectively, while the crosses on the lower left denote the median phot ometric error of eac h 
subsample. The extinction vector marks one magnitude of attenuation in the V band (ICardelh et al.lll989l) . 
Sources in the inner disk often diverge from the model track toward bluer values (top left), a likely effect of 
crowding. This could alternatively be attributed to differing filling factors in the gas, in which the emission 
originates (see Section 1575)) . Right: Two-dimensional color-density plot for the inner disk (blue), and outer 
galaxy (red). We omit areas that sample less than 15% of the 2D-histogram peak. The two distributions 
appear slightly offset, extending each other in opposite directions. The median colors (marked as stars) are 
separated by ~ 0.5 mag, possibly indicating a different age distribution for the two sub-populations in recent 
times (completeness affects the old end of the model track). 



occurs when binning LFs, while confirming the 
GielesI relation: the farther 'up' the function we 
move, the steeper the values become. 

To confirm this link to the iGieles we plot the 
slope of the LF as a function of brightness. We 
measure in 0.5 mag bins across the range shown 
in the LFs of Figure [T3] and find the slope to be- 
come steeper in proportion to the brightness. This 
evolution of the LF slope is shown in Figure 
Our results therefore suggest that all LFs might 
display substructure (also seen in the cumulative 
LFs of 
for the 



Bastian et al.ll2012l ). which would account 
Gieled relation, especially given its correla- 



tion with the range being fit. 

The LFs also provide a test of photometric com- 
pleteness. As outlined in Section [3741 we select in 



F336W but apply filters that expect the detec- 
tion of sources across the UBVI baseline. This 
implicitly includes the BVI bands in defining the 
completeness of our sample. In order to estimate 
this effect, we plotted the LFs by omitting the 
criteria pertaining to multi-band photometry and 
found the shape of each function to be unchanged 
by this process. The measured slopes were con- 
sistent within the errors with those listed in Ta- 
ble m We also examined the effect of potentially 
age-dependent reddening, by repeating the above 
exercise after limiting the age and extinction of the 
sample being fit. In all cases, the derived slopes 
were consistent with the fits of Table HI 
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Fig. 14. — Binned (left) and cumulative (right) luminosity functions (LF) of star cluster candidates in 
NGC 4041, arranged according to filter and Tgc (inner disk, outer galaxy). The lines (bands on the right) 
denote the range over which we fit the slope of each function, and we note slopes next to each pair of LFs. 
The U^se LF of the inner disk is shallower than other functions, which are generally consistent within the 
errors. A higher overall extinction would not explain this discrepancy. Slopes steeper than —2 are to be 
expected, given the bright absolute magnitudes being fit (Gieles 2010, see text for details). Binning smears 
the substructure evident in the cumulative functions and is subject to sampling errors. Right: LF slopes, 
measured on the cumulative functions of the cente r p a nel. Dashed lines indicate the 90% completeness 
limits. As expected by the literature review of Gidei (l2010l) . the slope becomes steeper with increasing 
brightness. We attribute this to the segmented structure of the function, as it gives rise to shallow slopes in 
the incomplete part, values consistent with the literature in the power-law segment, and very steep slopes 
once the upper truncation sets in. The increasing errors reflect the sampling (decreasing numbers). 



Table 4: Luminosity Function Fit Parameters and Outcomes. 



Filter Fit Range (in/out) Slope, unbinned (in/out) Slope, binned (in/out) 

H H 

F336W (U) [20.5,22.0] [21.5,23.5] 2.31 ± 0.24 2.56 ±0.15 2.35 ±0.10 2.56 ±0.10 

F450W (B) [20.5,22.0] [21.5,23.5] 2.94 ± 0.57 2.84 ± 0.56 3.21 ±0.12 2.94 ±0.22 

F606W (V) [20.5,22.0] [21.5,23.5] 3.13 ±0.40 3.06 ± 0.48 3.10 ±0.13 3.37 ±0.14 

F814W (I) [20.5,22.0] [21.5,23.5] 3.07 ± 0.37 2.93 ± 0.21 2.29 ±0.11 2.68 ± 0.05 

F336W (U) [20.0,23.6] [20.0,24.6] 2.01 ± 0.21 2.23 ± 0.23 1.60 ±0.04 1.84 ±0.06 

F450W (B) [20.0,23.6] [21.0,24.6] 2.46 ± 0.38 2.81 ± 0.16 2.23 ± 0.06 2.55 ± 0.04 

F606W (V) [20.0,23.6] [21.0,24.6] 2.37 ± 0.42 2.76 ± 0.25 1.90 ±0.07 2.60 ± 0.06 

F814W (I) [20.0,23.6] [21.0,24.6] 2.27 ± 0.43 2.57 ± 0.27 1.56 ±0.07 2.17 ±0.06 



Note. — The top and bottom tiers show fits to the 'smooth' and 'full' ranges (see text). The limits of the full range fits 
measure between 20 mag and the 90% completeness limits for the inner disk and outer galaxy (see Section [STTJ . 
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5.3. Mass Function 

In order to quantify the cluster mass distribu- 
tions in the inner and outer disk, we plot the Mass 
Function (MF) in Figure [15] as a set of cumulative 
distributions (counting from high-to-low mass). 
This plot will suffer from various incompleteness 
effects at different star cluster ages. On the young 
end we need to exclude short-lived, unbound stel- 
lar associations t hat masquerade as star clu sters at 
large distances ( Gieles fc Portegies ZwartI I2OIIV 
Toward older ages, the diagram is affected by 
the evolutionary fading and reddening of SSPs. 
This decreases the mass-to-light ratio, therefore 
increasing the mass required to detect a star clus- 
ter with increasing age. These evolving detection 
limits are represented by the dotted and dashed 
lines (inner/outer galaxy) on Figure[T5l and shape 
this parameter space. 

We therefore restrict the MF to cluster candi- 
dates aged between [10, 100] Myr. We note an 
offset of ~ 0.5 dex between the inner disk and 
outer galaxy, however the two functions display 
consistent slopes. Interestingly, the double-break 
structure reported for the LFs (Figure [13]) persists 
in the MF, which displays an incomplete part, a 
power-law segment, and an upper truncation. 

We investigate the origin of this segmentation 
by comparing to a set of power-law models, repre- 
sented by the various markings on Figure [TBI The 
slope is deduced as a = —2.3 through this pro- 
cess, consistent with the literature-standard index 
o f —2.0. In brief, w e follow the methodology of 



Bastian et al.l (|2012l ) and run a series of models 
of the MF with a sample size equal to the ob- 
served sub-populations, and assume a power-law 
distribution. The slope is derived after few trial- 
and-error iterations. The median of all models is 
shown as a solid line, and dashed and dotted lines 
indicate the range containing 50% and 90% of all 
modeled outcomes. The observed functions are 
seen to diverge from the power-law models at the 
low-mass end, ]bf and also the and high- mass end 
for the outer galaxy. While incompleteness will 
affect low-mass clusters, at the high-mass end we 
expect to detect the vast majority of sources. This 
divergence is therefore established as a physical ef- 
fect and the MF of the outer galaxy interpreted 
as a Schechter-type distribution with a trunca- 
tion at log(M,/Mo) w 4.4. This follows on the 
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Fig. 15. — Mass functions (MFs) of star clusters 
in the inner disk and outer regions of NGC 4041. 
The MFs cover only clusters aged 10 — 100 Myr 
and thus avoid fitting on young, unbound associa- 
tions, and incomplete samples at older ages. Both 
MFs show a composite structure in three parts: an 
incomplete low-mass section, followed by a power- 
law, which is truncated above log(MH,/M0) w 4.4 
in the outer galaxy. The lines represent a range 
of modeled pure power-law MFs of exponent —2.3 
(see text), which do not provide an adequate de- 
scription of the observed function for the outer 
galaxy, but describe the inner disk well. The light 
blue stars represent the three nuclear clusters in 
NGC 4041, as discussed in Section [4] 



analyses bv lGieled (|2009l) and|Larsen| (|2009( ) that 
favored Schechter MFs, and the truncation mass 
is in accord with th e value derived for M83 by 
Bastian et al.l (|2012[ ). As with the LF, the three- 
part structure of the MF suggests that binning 
acts to erase the information imprinted in the mass 
distribution. 

Interestingly, the MF of the inner disk is con- 
sistent with the power-law model, an effect poten- 
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tially related to the higher SFR in that region. 
The hght blue stars of Figure [TH top, represent 
the three nuclear region clusters discussed in Sec- 
tion H 

5.4. Age Distribution and Cluster Disrup- 
tion 

Combining mass with age we obtain the dia- 
gram of Figure ll6l The relative shift in mass noted 
above (Section I5.3[) is obvious here as a vertical 
shift according to different detection limits, while 
the age distributions appear largely similar and 
not unlike others studied in the past. The three in- 
ner disk outliers at high mass and young age repre- 
sent the double nuclear peak and another, nearby 
source. Assuming these are nuclear star clusters, 
we d o not expected them to follow norrnal sca ling 
laws (|Seth et al.l 120081 : IScott fc Grahamll2012h : as 
mixed stellar populations their mass-to- light ratios 
will be unlike those of SSPs. 

From the first age dex, log(r/?/r) G [6,7], we 
can make a rough assessment of the state of star 
formation in NGC 4041. While the outer galaxy 
might have recently produced a larger overall num- 
ber of clusters above the detection limit, the inner 
disk is consistently producing more massive clus- 
ters (the detection limit there is also brighter). 
The correlation between the mass of the m ost mas- 
sive cluster in a system and the SFR (jBastianI 
20081 ) indicates that the rate is higher in the in- 
ner disk, in accord with the value derived from the 
FUV fluxes (as noted in Section |44|) . 

Interpreting the distribution at older ages is 
more complicated. The formation of a star clus- 
ter does not guarantee its perpetual existence, as 
star clusters face many disrujjting processes over 
their early evolution (see iBastian fc Gieles .20081 
for an overview). In that way the presence of a 
cluster in a certain age bin is modulated not only 
by the fading lines discussed above, but also its 
survival as a physical entity. To decouple evo- 
lutionary fading from disruption, we can fit the 
age-mass distribution with a suite of models of 
mass-dependent cluster disruption. A dependence 
o n mass and environment was favored by the fits 
of lBastian et a l. (2012) and we follow this method- 
ology here. Other flav ors of cluster di ssolution - 
e. g. mass-independent |Fall et al.' '2005"), or purely 
environment-dependent (Elmcgrcen & Hunter 2010) 
- or combinations thereof, may also be valid. 
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Fig. 16. — Age-mass diagram for cluster candi- 
dates in the inner disk (blue circles) and outer 
regions (red triangles) of NGC 4041. The three 
young, massive sources in the inner disk (cyan 
stars) are located in the nuclear region, and may 
not follow the same scaling laws as SSPs. The 
lines show the 90% completeness limit translated 
into the mass of an Yggdrasil model across the 
given range of ages. This limit evolves accord- 
ing to the temporally increasing mass-to-light ra- 
tio of SSPs. The two sub-populations are offset 
by approximately 0.5 dex in mass, which follows 
on the relative shift of Figure [131 The age dis- 
tributions appear consistent and typical of their 
late-type host galaxy. 



We perform a maximum likelihood analysis of 
the age and mass distributions, shown in Fig- 
ur e [Tfl given th e clust er dis ruption formulations 
of iLamers et al.l (|2005l ) and iGieles] (|2009t ). This 
way we simultaneously compute and t^, the 
dissolution time of a 10^ cluster. We employ 
a Schechter function and the modeled evolution of 
an Yggdrasil SSP in Veoe? the limiting filter (see 
Section [37Tl Figure[3]). As with the LF analysis, we 
exclude sources younger than 10 Myr to avoid stel- 
lar associations, and take into account incomplete- 
ness at older ages (through the Yggdrasil model). 
This leaves the black dots of Figure [T71 which are 
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fit with models based on a range of variations of 
M, and ^4. 

The fits favor a t4 of a few hundred Myr in the 
outer galax y and a longer time in the inner disk, 
contrary to iLamers et al. ( 2005 ). who found en- 
vironment density to correlate inversely with t^. 
However, the small number of data-points in the 
inner disk provides broad confidence contours, ex- 
pressed as shallow topography, therefore no de- 
ductions should be made from this fit. The Af, 
seems to change significantly with local environ- 
ment (inner/outer galaxy), in a manner consistent 
with the M83 study from which this methodology 
is derived: denser regions display higher M*. 

5.5. Star Formation History 

A final test for a physical difference be- 
tween the inner/outer galaxy populations can 
be drawn from the cluster age distribution of 
Figure [18] (left panels, often referred to as a 
^dN/dt plot'), which counts the number of clus- 
ters surviving f rom every age-step of cluster 
formation (see iLamers et al.J _ l2005t iFall et al 



20051: IWhitmore et all 120071: IChandar et al.ll2010t 
Bastian et al. I2OIII and references therein) . When 
plotting all sources, as we do on the top panel, this 
diagnostic is affected by evolutionary fading and 
displays a power-law decline over time. Given a 
stable star formation history, setting a lower limit 
to the mass of a cluster sample will flatten the 
diagram to a certain age - the higher the mass 
cut, the farther back in time this diagnostic can 
reach. 

In the bottom panel of Figure [18] we perform 
this exercise for a mass cut that ensures com- 
pleteness to log(T/yr) = 8.5 (« 300 Myr). The 
precise mass of each cut (different for the inner 
and outer galaxy) is extracted from the age-mass 
diagram of Figure [H] as the intersection of the 
log(r/yr) = 8.5 line and the detection limit. The 
hashed regions are meant to guide the eye away 
from certain regions: the first two bins contain not 
only clusters, but also unbound associations; and 
the last bin is affected by incompleteness. That 
precaution leaves only three bins from which to 
make a deduction, hence we refrain from quanti- 
tatively characterizing the star formation history 
of NGC 4041 from this dN/dt plot. We note, 
however, that future SHUCS analyses focussed on 
nearby galaxies will not lack this diagnostic power. 



In previous Sections we have argued that bin- 
ning the typically small data-sets of extra-galactic 
star cluster populations erases information con- 
tained in the individual data-points. In order to 
investigate this effect in the age distribution, we 
add two more plots to the right panels of Fig- 
ure[l8] There we modulate the age distribution by 
the mass of each cluster, to obtain a mass output 
plot. The full populations show a slight dissimi- 
larity in their outputs over time, which, however, 
disappears when only plotting the complete sam- 
ple of high-mass clusters. 

5.6. A Possible Dynamical Event, Revis- 
ited 

In all, the study of star clusters in NGC 4041 
is in tune with the derivations of the traditional, 
broad-band diagnostics of Section [4] The color 
difference between inner and outer galaxy is mir- 
rored in the differing color, luminosity, mass, and 
age distributions of the corresponding star clus- 
ter populations. Unfortunately, the dN/dt plot is 
limited by the number of clusters available above 
the 90% completeness limit. Therefore, the past 
dynamical or accretion event deduced through 
multi-wavelength metrics above cannot be con- 
firmed through the cluster population. In follow- 
ing SHUCS investigations, however, we expect the 
smaller distances to the galaxies targeted to offer 
an opportunity to explore their dynamical histo- 
ries via the age distribution. 

6. Summary 

We have presented the observational and tech- 
nical setup of the Snapshot Hubble U-band Clus- 
ter Survey (SHUCS), comprising new F336W (U- 
band equivalent) observations of 10 galaxies with 
existing BVI coverage. We will use this new imag- 
ing to complete the UBVI baseline, as required 
for the accurate photometric age dating of star 
clusters. The survey is focussed primarily on star 
clusters in the first ^^Gyr of their evolution, and 
is aimed at understanding the formation and early 
evolution of these objects, as well as their utility 
as tracers of star formation history. In this paper 
we have also demonstrated the high scientific yield 
of combining star cluster statistics with broad- 
band, multi-wavelength metrics. We propose that, 
with large-scale information readily available from 
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Age [yr] [Mq] 

Fig. 17. — Disruption time for a 10^ M© cluster (^4) versus Schechter Mass (M*, where the truncation sets 
in) for the inner disk and outer regions of NGC 4041. This plot uses the age-mass diagram of Figure [TBI 
and completeness (based on an Yggdrasil SSP model) to estimate the two quantities through a series of 
maximum- likelihood tests. Similarly to the MFs of Figure [121 only the black dots are used, in order to avoid 
incompleteness effects and short-lived stellar associations. We find values consistent with the literature in 
the outer galaxy, while the atypically high Af* derived for the inner disk might be affected by crowding and 
incompleteness. This is evident from the shallow gradient of the contours, i. e. the large range of values that 
fit the data well. 



GALEX, SDSS, 2MASS, IRAS, and other surveys, 
this two-pronged approach is not only advanta- 
geous, but imperative in achieving a well-rounded 
understanding of each cluster population at hand. 

This combined approach was demonstrated 
through the analysis of the large-scale properties 
and stellar populations of NGC 4041, a massive 
galaxy at the upper tier of the SHUCS distance 
scale {d > 15 Mpc). This system features a com- 



plex physical and dynamical structure, expressed 
as a discontinuous brightness profile perhaps rem- 
iniscent of M64 from a different vantage point. We 
find a strong distinction in the colors and masses of 
star clusters when applying a cut in galactocentric 
distance where this break occurs, reminiscent of 
the color and age distributions in interacting and 
merging systems. Most notably, we find the inner 
disk to have been forming more massive clusters 
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Fig. 18. — The cluster age distribution histograms {'dN/dV plots) of the left panels count the number 
of surviving clusters per time interval over the past Gyr. Isotropic axes show an x = —y slope in the 
magnitude-limited sample on the top, as expected from the literature. Applying a mass cut should flatten 
the dN/dt of a galaxy with a stable star formation history. In this case, the slope lessens, but is difficult 
to interpret given the small number of data-points. In this work we have argued that binning might erase 
information, therefore we also present unbinned data in the cumulative age distribution on the (right). The 
two sets of plots show consistent results. 



(M > 10** M0) in the past w 100 Myr, despite 
the overall redder appearance of its subpopula- 
tion, which indicates an older mean age. This 
is in tune with our analysis of archival GALEX 
FUV imaging, which reveals a higher SFR in 
the inner disk than the outer galaxy, similar to 
post-starburst systems. This could be linked to 
the theoretical expectation of a higher star for- 
mation efficiency and a higher fraction of stars 



form ing in clusters i n regions of higher gas den- 
sity (|Kruiissenll2012l ). We also discovered a tidal 
feature, which, combined with information in the 
literature, strongly favors an accretion event in 
the recent past 100 Myr) as the origin of the 
segmented morphology of NGC 4041. 

Throughout the galaxy we found a truncation 
in the star cluster mass function, in accord with 



the recent studies ( Larsen 120091: iGielesI [Sool 
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This truncation occurs at a higher mass in the 
inner disk, as expe cted by these recent resuhs 
( Bastian et aLllioiih . Furthermore, we found both 
the mass and luminosity function to break down 
into three segments when binning is abandoned 
in favor of cumulative distributions: one dom- 
inated by incompleteness, one following the fa- 
miliar power-law shape, and one encompassing a 
truncation at the upper end. Our results strongly 
advocate for the use of all available information, 
rather than binning, when characterizing the age, 
mass, and luminosity distributions of star cluster 
populations. In this work we have sought that re- 
sult through the use of cumulative distributions. 

The strength of the survey-at-large derives 
from the availability of deep 1/33^ observations at 
the highest resolution available (currently HST- 
WFC3), enabling the precision age dating of hun- 
dreds of clusters per galaxy - potentially thou- 
sands in nearby galaxies. The full survey will 
consist of the analysis of the cluster populations 
of ten late-type galaxies of various morphological 
and spectroscopic types. In addition to the un- 
precedented statistical value of this analysis, we 
expect to focus a few works on individual sys- 
tems, namely the ongoin g merger in NG C 2146 
(see the pilot study of Adamo et al.ll2012l ). grand 
design spiral NGC 2997, and NGC 247, where we 
will contrast the star formation history derived 
from field stars to that derived from star clusters. 
Combined, the individual studies will help: 

1. Search for a characteristic value (M*) in the 
cluster mass function. 

2. Empirically constrain star cluster disruption 
laws. 

3. Constrain the fraction of stars that form in 
clusters. 

4. Plot the star/cluster formation histories of 
SHUCS galaxies over the past Gyr. 

5. Determine whether cluster size is affected by 
local conditions. 

6. Search for environmental dependencies in all 
of the above. 
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